Pattern Formation and a Clustering Transition in Power-Law Sequential Adsorption 
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A new model that describes adsorption and clustering of par- 
ticles on a surface is introduced. A clustering transition is 
found which separates between a phase of weakly correlated 
particle distributions and a phase of strongly correlated dis- 
tributions in which the particles form localized fractal clus- 
ters. The order parameter of the transition is identified and 
the fractal nature of both phases is examined. The model 
is relevant to a large class of clustering phenomena such as 
aggregation and growth on surfaces, population distribution 
in cities, plant and bacterial colonies as well as gravitational 
clustering. 



Many of the growth and pattern formation phenom- 
ena in nature occur via adsorption and clustering of par- 
ticles on surfaces [0j2). The richness of these phenom- 
ena may be attributed to the great variety of structures 
and symmetries of the adsorbed particles and substrates. 
Nonequilibrium growth models often give rise to frac- 
tal structures, which are statistically self similar over a 
range of length-scales |^ . In a large class of surface ad- 
sorption systems, the dominant dynamical process is the 
diffusion of the adsorbed particles which hop randomly 
on the surface until they nucleate into immobile clusters 

. The formation of fractal clusters, which are typical 
in these systems can be described by the diffusion lim- 
ited aggregation (DLA) process Q . In DL A a cluster of 
particles grows due to a slow flux of particles which dif- 
fuse as random walkers until they attach to the cluster. 
The model describes a great variety of aggregation pro- 
cesses such as island growth in molecular beam epitaxy 
electrodeposition, viscous fingering, dielectric break- 
down, and various biological systems In many other 
physical systems, once an adsorbed particle sticks to the 
surface it becomes immobile. These systems can be de- 
scribed by random sequential adsorption (RSA) processes 
Within the RSA processes, one should distinguish 
between systems in which particles cannot overlap and 
systems in which they can adsorb on top of each other. 
Systems in which particles cannot overlap tend to reach 
a jamming limit, in which the sticking probability of a 
new particle vanishes . Models which allow multilayer 
growth describe a large class of physical systems, includ- 
ing deposition of colloids, hquid crystals polymers 
and fiber particles [p|Jlo[| . Recently, the case of power law 
distribution of particle sizes was studied both for uncor- 
related adsorption fll]] and for non-overlapping particles 



|]l2| . In the case of uncorrelated adsorption it was found 
that the boundary of the particle clusters is fractal . 
For non-overlapping particles, it was found that the area 
that remains exposed is fractal [ p^ . 

Models that describe growth dynamics have been em- 
ployed in recent years in a vast range of scientific fields as 
diverse as city organization and growth city and 

highway traffic ||l^ and growth of bacterial colonies . 
A common feature is the tendency of the basic objects to 
form clusters of high density (typically of fractal shape), 
surrounded by low density areas or voids. Other exam- 
ples of clustering appear in the distribution of mass in the 
universe , in dissipative gases and granular fiow as 
well as in step bunching on crystal surfaces during growth 
p9| and due to electromigration |2^. The phenomenon 
of cluster formation is therefore generic in a broad class 
of systems in spite of the fact that the pattern forming 
dynamical processes may vary substantially from system 
to system. This richness of clustering phenomena is not 
yet fully backed-up by appropriate models. 

In this Letter we introduce a new model: the Power- 
Law Sequential Adsorption (PLSA) model which de- 
scribes a variety of surface adsorption and clustering pro- 
cesses. This model leads to a rich variety of structures, 
many of which are fractal, which mimic the experimental 
morphologies found in the examples cited above. In par- 
ticular, it exhibits a unique clustering transition which 
separates between a weakly correlated phase in which 
the adsorbed particles are distributed homogeneously on 
the surface and a strongly correlated phase in which they 
form clusters. 

In the PLSA model circular particles of diameter d are 
randomly deposited on a two dimensional (2D) surface 
one at a time. The deposition process starts from an ini- 
tial state where there is one seed particle on the surface. 
The sticking probability < p < 1 of a newly deposited 
particle is determined by the distance r from its center 
to the center of the nearest particle which is already on 
the surface. This probability is given by 



p(r) 



1 r <d 

{d/rY r>d, 



(1) 



where the exponent a > is a parameter of the model. 
The model thus exhibits a positive feedback clustering, 
like many of the clustering phenomena listed above. The 
random deposition process is repeated until the desired 
number of particles, M, stick to the surface. Since the 
sticking probability is given by a power law function, 
which is of a long range nature, this process should have 
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been studied, in principle, in the infinite system limit. 
Numerical simulations, however, are done on a finite sys- 
tem of area Lx L, where L/d:^ 1. The coverage is given 
by 77 = A/L'^ where A is the total area covered by the 
particles. Also, let 770 = 7r(d/2)2 • M/L^. 

The limit of uncorrelated adsorption, in which the 
sticking probability is p = 1 uniformly, is obtained for 
a = Q. This limit was studied recently using fractal anal- 
ysis, and the box-counting and Minkowski functions were 
calculated analytically |21|] . It was found that for a range 
of low coverages, apparent fractal behavior is observed 
between physically relevant cutoffs. The lower cutoff rp 
is given by the particle diameter = d while the upper 
cutoff ri is given by the average gap between adjacent 
particles, namely ri = p^^/^ — d, where p = MjL?' is the 
particle density. 

The limit of strongly correlated adsorption is obtained 
for a — > 00. In this case only a single, connected cluster 
is generated on the surface. The perimeter of this clus- 
ter grows slowly when new particles are deposited on its 
edge, while it becomes more dense inside p|,p^. 

We will now examine the morphological properties of 
the configurations of adsorbed particles for the full range 
of < a < cx) using fractal analysis. For this analysis we 
use the box-counting (BC) procedure in which one cov- 
ers the plane by non-overlapping boxes of linear size r. 
The box-counting function N{r) is obtained by counting 
the number of boxes which have a non-empty intersec- 
tion with the (fractal) set. A fractal dimension Dec, is 
declared to prevail at a certain range of length-scales, if 
a relation of the type NBc{f) ^ r'^^^c holds, or equiv- 
alently, if the slope of the log-log plot 

(2) 



Dbc = -slope {logr, log[iVBc(f)]} 



is found to be constant over that range. 

Two configurations of particles, randomly deposited 
and adsorbed according to Eq. (Q) onto the unit square 
(L = 1), are shown in Fig. 1 for 770 — 0.01. For a — 1.5 
the particle distribution exhibits local density fluctua- 
tions but on larger scales it is rather homogeneous and 
extends over the entire system [Fig. 1(a)]. For a = 2.5 we 
observe a strongly clustered structure [Fig. 1(b)]. This 
structure resembles the set of turning points of a Levy 
flight random walker [^^. In fact, a Levy flight corre- 
sponds to the special case in which the sticking proba- 
bility of the next deposited particle depends only on the 
position of the latest particle adsorbed on the surface. 
Unlike Levy flights which typically describe dynamic be- 
havior our model describes clustering in spatial struc- 
tures. It is also related to other models of spatial struc- 
tures such as the continuous percolation model which is 
approached when the interaction is suppressed, at a — > 0. 
Another related model, which describes the growth of a 
percolation cluster and exhibits power-law correlations 
between growth sites is presented in Ref . . 

The box-counting functions for the conflgurations gen- 
erated by the PLSA model are shown in Fig. 2. It is ob- 
served that for a < 2 the box-counting function resembles 



the shape obtained for the uncorrelated case [£l[. This 
indicates that the basic features of the model studied in 
Ref. [H^ are maintained not only for short range corre- 
lations but for an entire class of long range correlations. 
The box-counting function for a > 2 exhibits a nearly 
linear behavior for the entire range from the particle size 
to the cluster size. 

To obtain the fractal dimensions of the sets from the 
box-counting functions one should identify the relevant 
range of length-scales over which the linear fit should be 
done. For the weakly correlated distributions the relevant 
range of length scales spans between ro = d and ri = 



-1/2 _ d 



For the strongly correlated distributions 



P _ 

where clusters are formed, the relevant range is limited 
from above by the linear size of the entire cluster. The 
quality of the linear fit is measured by the coefficient of 
determination TZ^ [^. In both cases, given a desired 
value of 7^^ one can further narrow the range within the 
cutoffs described above to find the broadest range (ro, ri) 
within which the linear regression maintains the given 
value of 7^2 ||. 

The fractal dimension Z? as a function of a is shown 
in Fig. 3 for rjo = 0.1, 0.01 and 0.001. Two domains 
are identified: a plateau of low dimension for the weakly 
correlated case and a plateau of high dimension for the 
strongly correlated case. 

Consider a seed particle located at the origin of an 
infinite plane. Particles are randomly deposited one at 
a time according to the PLSA model until one particle 
sticks to the surface. Consider the probability that the 
distance r between the first particle which sticks and the 
seed particle at the origin will be larger than some value 
rf where rf > d. This probability is given by: 



P(r > rf) = 



l°^27rr{lrdr 
f^2Trrdr + f^2Trr{!i)"dr' 



(3) 



One readily verifies that for a < 2 the probability P{r > 
rf) ^ I for any finite rf. For a > 2, on the other hand, 
this probability is given by 



P{r >rf)^^ f A 
a \rf 



(4) 



Therefore, in the infinite system limit of the weakly cor- 
related phase {a < ac) the probability that the next 
particle will stick within any finite distance from an ex- 
isting cluster is zero [^ . In the strongly correlated phase 
(a > ac), the probability that the next particle will stick 
within a finite distance rf from the cluster can be made 
arbitrarily close to one, by an adjustment of accord- 
ing to Eq. (Q) [^^. In general, the value of ac for which 
the clustering transition takes place is equal to the space 
dimension. 

The order parameter of the clustering transition is 1^ = 
{rjQ — 77)7770, namely the fraction of the total area of the 
adsorbed particles which is lost due to overlap. Consider 
a finite number M of particles of diameter d — 1 adsorbed 
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on the surface in the infinite system Umit (L oo). For 
a < 2, in this low coverage hmit overlaps are negligible 
and y = 0. For a > 2 clusters become more dense and 
overlaps more dominant as a increases. Our numerical 
studies are done on a finite system of size L = 1 for a 
range of coverages. The order parameter as a function 
of a, for 7^0 — 0.01, is shown in Fig. 4. 

To examine the critical behavior in the infinite system 
limit we performed analytical calculations in one dimen- 
sion (ID). In ID the configuration is fully specified by 
the ordered list of M — 1 gaps between the M parti- 
cles. For the ID case we have obtained the critical ex- 
ponent (3 for the order parameter V ^ {a — ac)^ in the 
L —>■ oo limit by constructing the probability distribu- 
tion Pi, i — 0, . . . ,M that the next particle that sticks 
will stick within the gap gi (where go and gM are the two 
semi- infinite gaps on both sides). The overlap was then 
calculated as a weighted average over all gaps. The re- 
sult we obtain is that the exponent /3 = 1. We also found 
that the fractal dimension exhibits critical behavior of 
the form D ~ (a — ttc)''', where ac = 1 with the exponent 
7=1. For the weakly correlated phase at a < 1 the 
fractal dimension in the L ^ oo is -D = while for a > 2 
the dimension is D = 2. 

In summary, we introduced a new model for random 
sequential adsorption, characterized by a power law dis- 
tribution of sticking probabilities. This model exhibits 
a continuous phase transition between weakly correlated 
adsorption, in which the particle distribution is homoge- 
neous on large scales and extends over the entire system, 
and strongly correlated adsorption in which a highly clus- 
tered structure is generated. We thus identified a broad 
class of distributions which maintain the basic proper- 
ties of the weakly correlated random structures studied 
in Ref. |^ and found the borderline between this class 
and the class of strongly correlated structures which ex- 
hibit clustering phenomena. The model should be useful 
in the study of a great variety of clustering problems. 
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FIG. 1. Particles adsorbed on the surface of a unit square {L = 1) according to the PLSA model for rjo = 0.01: (a) for a = 1.5 
we observe density fluctuations at small scales, however at larger scales the distribution is rather homogeneous and extends 
over the entire system; (b) for a — 2.5 wc observe a strongly clustered structure and vacant area elsewhere. The number of 
particles in both (a) and (b) is 3184 and their diameter is d = 0.002. 

FIG. 2. The box-counting function for four configurations with rjo = 0.01, and a = (empty circles), 1.5 (full circles), 2.5 
(empty squares) and 3.5 (full squares). It is observed that for a < 2 the shape of this function resembles that of the uncorrclatcd 
case. For a > 2, where strongly clustered distributions arise, there is a broad scaling range. The units are dimensionless and 
the logarithms are in base 10. 

FIG. 3. The fractal dimension of the configurations produced by the PLSA model as a function of a (in dimensionless units) 
for r]o = 0.1 (empty circles), 0.01 (full circles) and 0.001 (empty squares). A plateau of low dimension is found for the weakly 
correlated limit and a plateau of high dimension for the strongly correlated limit. 

FIG. 4. The order parameter of the clustering transition, V = (rjo — ri)/rio, is shown as a function of a for rjo = 0.01. It 
represents the fraction of the total area of the adsorbed particles which is lost due to overlap. For r)o <^ 1, this order parameter 
vanishes for a < 2 and increases above a = 2. 
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